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ABSTRACT 
The  cold  plasma  conductivity  tensor  is  employed  to  explore  the 
problem  of  wave  propagation  in  a  plasma  withovit  the  use  of  the 
geometrical  optics  approximation  or  a  plane  stratified  model.  After 
the  introduction  of  the  appropriate  electromagnetic  potentials,  the 
wave  propagation  problem  reduces  to  a  sixth  order  self-adjoint  set  of 
partial  differential  equations  derivable  from  a  self-adjoint 
variational  principle.  The  usual  lower  and  upper  hvhrid  resonances  as 
well  as  the  plasma  resonances  can  he  identified  as  surfaces  across 
which  the  svstem  of  equations  changes  from  elliptic  to  partially 
elliptic  and  partially  hyperbolic.  In  addition,  the  ion  and  electron 
cyclotron  resonance  surfaces,  near  which  the  cold  plasma  approximation 
breaks  down,  are  also  transition  surfaces. 

When  the  underlying  plasma  equilibrium  is  axisvmmet r i c ,  although 
the  waves  need  not  be  symmetric,  a  simpler  set  of  potentials  is  given 
for  which  the  wave  propagation  problem  reduces  to  a  fovirth  order 
system.  The  equations  are  again  derivable  from  self-adjoint 
variational  principles.  The  bulk  of  the  paper  considers  the 
axisymmetric  equilibrium  case.  For  waves  prooa^ating  at  about  the  ion 
cyclotron   frequency,   it   is   proposed   that   the   electric   field   is 
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continuous  and  Its  normal  derivative  is  discontinuous  across  the  ion 
cyclotron  resonance  surface.  With  this  stipulation  wave  patterns  and 
heating  rates  may  be  calculated  in  this  freqviency  range.  It  is 
indicated  how  non-standard  such  a  problem  is.  At  frequencies  on  the 
order  of  the  ion  plasma  frequency  and  including  the  lower  hybrid 
resonance  frequency  a  different  situation  pertains.  First,  a  different 
geometric  optics  model  is  proposed  in  which  certain  spatial  gradients 
of  equilibrium  quantities  appear.  The  origin  of  the  difference  is  the 
connection  of  the  value  of  the  small  geometrical  optics  expansion 
parameter  to  the  magnitude  of  various  elements  of  conductivity  tensor. 
The  usual  expansions  assume  that  all  the  elements  of  the  tensor  are  of 
order  one  in  the  geometrical  optics  parameter.  Next,  the  equations  are 
studied  near  the  resonance  surface.  It  is  argued  that  if  the  resonance 
surfaces  are  not  the  same  as  flux  surfaces  then  heating  occurs  only  at 
the  points  of  tangency  of  flux  surfaces  and  the  resonance  surface.  In 
this  case  the  system  is  compared  with  the  Tricomi  equation  of  transonic 
flow.  If  the  resonance  surfaces  and  flux  surfaces  coincide  then 
heating  is  possible  on  the  entire  resonance  —  or  flux  —  surface.  The 
great  complexity  of  the  general  wave  propagation  and  heating  problem  in 
this  frequency  range  is  described. 
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1.  INTRODUCTION 

The   study  of  wave   propagation   originated  decades  ago  with  the 

examination  of  radio  wave  propagation  in  the  ionosphere.   Once  the  cold 

plasma  conductivity  tensor  is  developed,  the  problem  reduces  to  a  study 

of  Maxwell's  equations  with  the  specified  conductivity  tensor.   Ry  far, 

the   bulk   of   the  work  done  in  this  field  has  been  based  on  geometric 

optics  or  plane  stratified  medium  approximations.   Neither   of   these 

approximations   is  particularly  valid  in  many  typical  laboratory  plasma 

applications   at   moderately   low   frequencies   and   moderately   long 

wavelengths.    Additionally,   the   appearance   of   "resonances"  and  the 

associated  plasma  heating  brings  into  question   the   relevance   of   the 

simple   cold  plasma  model.   Here,  we  intend  to  examine  fairly  generally 

the   consequences   of   the   cold   plasma   model   for   wave   propagation 

problems.   We   intend   to   compare   the  situation  in  a  plasma  with  the 

better  known  problems  of  scattering  of  waves  in  electromagnetic   theory 

or   in  quantum  mechanics.    Our  attack  on  this  problem  is  based  on  an 

examination  of   the   system  of   partial   differential   equations   that 

defines   the  wave   propagation   problem.    We  will  not  be  able  to  give 

definitive  answers  to  many  of  the  questions  we  raise,  as   the   problems 

appear   far   beyond   current  knowledge.    We   will   show  that  waves  in 

plasmas  behave,  in  many  cases,  quite  unlike  the  much  simpler  version  of 

waves   in  ordinary  electromagnetic  theory  or  in  quantum  mechanics.   We 

believe  that  the  questions  we   raise   must   be   addressed   directly   to 

obtain  reliable,  quantitative  measures  of  ion  cyclotron  or  lower  hybrid 

wave  heating  and   model   patterns   in   a   plasma   device.    \nv      direct 

numerical  calculations  face  immediately  the  problems  posed  here. 

We   employ  Maxwell's   equations  with  the  cold  plasma  conductivity 
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tensor  and  we  assume  time  harmonic   dependence   of   H   and   B.    As   is 

well-known   this   set   of   partial   differential   equations   is   highly 

degenerate  since  the  operator  Vx   has   the   large   null   space   of   all 

gradients.    Thus,  the  usual  prescriptions  of  how  to  determine  the  tvpe 

of  a   system,   elliptic  hyperbolic,   parabolic,   or   mixed,   gives   no 

Information.   We   are  unable  even  to  give  the  order  of  the  system.   We 

show  that  the  introduction  of  the  usual  electromagnetic   potentials   in 

an  appropriate   gauge   generates   a   mathemat icallv  sensible  system  of 

equations,  in  fact,  the  system  is  self-adioint  and  mav  be  derived   from 

a   variational   principle.    The   system,   not   surprisingly,   can   be 

identified  as  being  sixth  order.   We  can  now  determine  the  tvoe  of   the 

system.    Typicallv,   time   harmonic  wave   prooagation   or   scattering 

problems  are  represented  by  elliptic  systems.   We   find   that   for   the 

cold   plasma  this  is  not  necessarily  the  case.   We  find  that  the  system 

always  has  some  elliptic  characteristics,  but   there   are   additionally 

transitions    from   elliptic   to   hyperbolic   behavior   at   the   plasma 

frequency  resonance  surface,  lower  and  upper  hvbrid  resonance  surfaces, 

and   at   the  fundamental  ion  and  electron  cvclotron  resonance  surfaces. 

The  usual  introduction  of  these  resonances  depends  on  properties  of   an 

infinite  homogenous  medium  or  possiblv  a  stratified  medium.   We  now  see 

that  the  surfaces  are  related  to  an  essential  pronertv  nf   the   system. 

The   nature   of   the   solution   to   wave  prooagation  problems  near  such 

surfaces  is  one  of  the  major  topics  of  this  work.    We   note   that   the 

cyclotron   resonances   are  fundamentally  different  from  the  others.  \t 

cyclotron   resonance   elements   of   the   condnctivitv    tensor   become 

unbounded,   clearlv   indicating   a  failure  of  the  cold  plasma  model,  at 

the  other  resonances  there  is  no  obvious   failure   of   the   model.    We 
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shall  see  that  the  solutions  are  fundamentally  different  near  cyclotron 
resonance  and   lower  hybrid   resonance,   as  well.    We   believe   that 
cyclotron   resonance  and  lower  hybrid  resonance  refer  to  fundamentally 
different  phenomena. 

The  preceding  results  were  obtained  for  the  fully  general, 
three-dimensional  case.  In  many  cases  one  is  interested  in  wave 
propagation  in  a  system  in  which  the  underlying  equilibrium  is 
axisymmetric  although  the  waves  lack  this  symmetry.  For  this  situation 
we  can  give  a  better  formulation  of  the  wave  propagation  problem.  We 
give  two  different  representations,  depending  on  whether  the  wave 
fields  preserve  the  axisyraraetry  of  the  equilibrium  or  they  destroy  it. 
In  both  cases  the  electromagnetic  fields  are  represented  in  terms  of 
two  potentials.  The  equations  are  of  fourth  order  rather  than  sixth 
order  as  in  the  general  three-dimensional  case.  We  are  unsure  whether 
the  change  in  order  comes  from  a  poor  formulation  in  the 
three-dimensional  case  or  partial  integration  of  the  system  in  the 
axisymmetric  case.  The  characteristic  surfaces  for  both  of  our  fourth 
order  systems  are  the  same  and  correspond  to  two  elliptic  (imaginary) 
characteristics  plus  two  others  that  may  be  either  elliptic  or 
hyperbolic.  The  surfaces  of  transition  of  type  are  not  the  same  as  in 
the  three-dimensional  case,  although  the  upper  and  lower  hybrid 
resonance  transitions  are  unaffected.  The  systems  of  partial 
differential  equations  are  derivable  from  self-adjoint  variational 
principles. 

The  elements  of  the  cold  plasma  conductivity  tensor  are  very 
different  in  magnitude.  In  order  to  proceed  further  we  must  select  a 
frequency  regime  and  other  physical  parameters.   For  laboratory  plasmas 
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we  assume  that  the  electron  cyclotron  and  plasma  frequencies  are 
comparable  and  we  examine  separately  waves  at  about  the  ion  cyclotron 
frequency  and  the  ion  plasma  frequency.  The  underlying  small  parameter 
to  which  we  compare  others  is  the  electron-ion  mass  ratio.  We  consider 
first  waves  at  about  the  ion  cvclotron  frequency.  The  situations  in  a 
mirror  machine  and  a  Tokamak  are  somewhat  different  and  it  is  easier  to 
describe  them  separately.  In  a  mirror  machine  there  is  an 
elliptic-hyperbolic  transition  at  the  fundamental  ion  cyclotron 
resonance  surface.  The  system  is  elliptic  for  frequencies  above  the 
ion  cyclotron  resonance  frequency  and  hyperbolic  below.  Near  the 
resonance  surface  the  elements  of  the  conductivitv  tensor  are  unbounded 
and  the  cold  plasma  model  fails.  In  another  paper  (1)  we  have  analyzed 
this  situation  with  the  help  of  a  better  model  without  singularities, 
the  Vlasov  equation,  and  we  have  shown  that  the  solution  to  the  wave 
propagation  problem  can  be  described  in  terms  of  the  cold  plasma  model 
with  the  additional  hypotheses  that  the  electric  field  is  continuous 
across  the  resonance  surface,  but  that  its  normal  derivative  suffers  a 
discont inuitv.  Additionally  one  can  calculate  the  heating  in  terms  of 
a  simple  surface  integral.  Thus,  the  solution  to  the  wave  propagation 
and  energy  absorption  problem  becomes  the  calculation  of  the 
electromagnetic  fields  for  a  system  which  in  some  regions  is  elliptic 
and  in  other  regions  is  partly  elliptic  and  partly  hyperbolic.  Across 
the  transition  surface  the  fields  have  specified  discontinuities  in 
their  normal  derivatives.  This  problem  is  quite  non-standard,  not  only 
in  the  mixed  type  of  the  system,  but  also  in  the  appearance  of  jump 
conditions.    Presumably   the   self-adjoint  variational  principles  will 
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allow  direct  calculation  of  the  fields,  but   it   should  be   emphasized 
that  the  problem  is  far  from  the  ordinary. 

In  a  Tokamak  the  ion  cyclotron  resonance  problem  is  somewhat 
different.  Formally,  the  system  of  equations  exhibits  no  singularity 
at  the  ion  cyclotron  resonance  although  there  is  an  elliptic-hyperbolic 
transition  very  close  to  the  ion  cyclotron  resonance  surface.  The 
elliptic-hyperbolic  transition  occurs  within  the  region  in  which  the 
cold  plasma  model  likely  fails.  Although  we  have  not  carried  through 
the  detailed  analysis  we  believe  that  analysis  similar  to  what  has  been 
done  for  the  mirror  machine  would  gives  similar  results.  We  expect  the 
fields  to  be  described  by  the  cold  plasma  model  except  for  jump 
conditions  across  the  transition  surface.  It  is  interesting  to  note 
the  appearance  of  the  result  that  the  component  of  the  time  varying 
electric  field  parallel  to  the  steady  magnetic  field  is  small.  In  our 
treatment  the  result  follows  from  the  representation  of  the  fields  in 
terms  of  potentials  and  is  not  a  property  of  the  detailed  solutions  for 
the  potentials.  This  result  holds  for  both  cases  of  mirrors  and 
Tokamaks. 

When  we  turn  to  the  study  of  waves  propagating  at  frequencies  on 
the  order  of  the  ion  plasma  frequency  and  waves  at  lower  hybrid 
frequencies  we  find  a  much  more  complex  problem  with  many  open 
questions.  We  consider  in  detail  only  axisymmetric  modes  in  a  Tokamak 
although  the  results  holds  for  non-axisymmetric  modes  in  a  Tokamak  as 
well.  We  have  not  explored  the  corresponding  problems  in  a  mirror. 
The  elliptic-hyperbolic  transition  occurs  exactly  on  the  lower  hybrid 
resonance  surface  and  there  are  no  other  potential  singularities  of  the 
system  in  this  frequency  range.   Before  we  address  the   solutions   near 
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the  resonance  surface  we  consider  the  usual  geometric  optics  model 
approximation  for  the  system.  A  direct  examination  of  our  system  of 
equations  shows  that  we  obtain  in  the  geometric  optics  limit  a 
different  set  of  equations  from  those  one  would  find  by  the  replacement 
of  V  by  ik  in  Maxwell's  equations.  Specifically,  we  find  that  spatial 
gradients  of  equilibrium  quantities  appear  in  our  svstem.  We  found  a 
similar  phenomenon  in  a  study  with  another  author^-^  of  lower  hybrid 
waves  in  a  perpendicularly  stratified  medium.  The  origin  of  the 
difference  is  quite  clear.  In  a  physically  reasonable  modelling  of  a 
geometrical  optics  approximation  to  lower  hybrid  waves, the  geometrical 
optics  expansion  parameter  appears  explicitly  and  non-trivially  in  the 
conductivity  tensor.  The  usual  geometrical  optics  treatments  omit  such 
dependence.  Thus,  even  where  geometrical  optics  is  a  good 
aDproximation  we  expect  significantly  different  oropagation 
characteristics  from  those  usually  given. 

Finally  we  consider  wave  propagation  near  the  lower  hybrid 
resonance  layer.  We  show  first  that  the  governing  differential 
equations  are  very  different  depending  on  whether  or  not  the  plasma  is 
perpendicularly  stratified.  If  the  medium  is  not  peroendicularlv 
stratified,  the  equations  in  normal  form  are  essentially  the  Tricomi 
equation  of  transonic  flow,  where  the  resonance  surface  is  iust  the 
sonic  transition  surface,  and  Laplace's  equation.  It  has  been 
well-established  that  the  Tricomi  equation  with  appropriate  boundary 
data  has  smooth,  well-behaved  solutions . ^ ^ ^  Thus,  we  expect  the  wave 
propagation  problem  to  have  smooth  solutions  without  any  energy 
absorption  across  most  of  the  resonance  surface.  Near  those  points  at 
which  the  resonance  surface  is  tangent  to  a  flux  surface,  the  system  is 
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unlike  the  Tricomi  equation  and  our  argumentation  fails,  and  we  believe 
that  energy  absorption  is  possible.  Similarity  solutions  for  a  model 
partial  differential  equation  with  the  properties  of  lower  hybrid  waves 
show  just  these  phenomena  of  no  singularities  or  heating  on  flux 
surfaces  or  resonance  surface  except  where  they  are  parallel.  When  the 
plasma  is  perpendicularly  stratified,  the  situation  is  quite  different 
and  we  expect  heating  along  the  entire  resonance  —  or  flux  —  surface 
exactly  as  in  our  earlier  work.^'^^  We  cannot  guarantee  the  surety  of 
these  results,  but  we  believe  thev  correspond  to  what  is  generally 
expected  of  such  systems  of  equations.  Even  accepting  the  validity  of 
our  assertions  we  are  not  yet  ready  to  formulate  the  energy  absorption 
model  near  those  points  at  which  flux  surfaces  are  tangent  to  resonance 
surfaces.  The  full  solution  of  a  wave  propagation  and  energy 
absorption  problem  requires  the  resolution  of  this  matter. 

The  organization  of  this  paper  follows  the  subject  matter  of  this 
introduction.  The  next  section  shows  how  to  formulate  Maxwell's 
equations  with  the  cold  plasma  conductivity  tensor  in  a  sufficiently 
explicit  manner  as  to  allow  the  determination  of  the  mathematical  type 
of  the  set  of  oartial  differential  equations.  The  surfaces  on  which 
the  type  changes  are  found  and  shown  to  be  the  usual  cvclotron,  hybrid 
and  plasma  resonance  surfaces.  The  third  section  shows  how  to 
introduce  a  simpler  electromagnetic  field  representation  in  terms  of 
two  potentials  for  the  case  that  the  underlying  plasma  equilbrium  is 
axisymmetric.  The  mathematical  type  of  the  system  is  also  given  as 
well  as  the  surfaces  of  elliptic-hyperbolic  transition.  In  preparation 
for  subsequent  material  the  scaling  of  the  various  small  parameters  in 
the  electron  to  ion  mass  ratio  is  given.   In  the   Fourth   section   wave 
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propagatlon  at  about  the  ion  cyclotron  frequency  is  studied.  The 
elliptic-hyperbolic  transition  is  essentially  the  ion  cyclotron 
resonance  and  the  behavior  near  the  resonance  surface  is  described. 
The  last  section  treats  waves  at  about  the  ion  plasma  frequency.  A 
modified  geometrical  optics  which  includes  some  spatial  gradient 
effects  is  proposed.  It  is  argued  that  plasma  heating  occurs  only  near 
those  points  at  which  the  resonance  surfaces  are  parallel  to  the  flux 
surfaces. 
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II.   SOME  GENERAL  PROPERTIES  OF  THE  COLD  PLASMA  MODEL 

In  this  section  we  give  some  standard  material  on  Maxwell's 
equations  in  a  dielectric  medium  and  we  explore  some  properties  of  the 
system.  Our  goal  is  to  uncover  enough  Information  about  the  equations 
so  that  we  may  propose  plausible  methods  for  determining  solutions 
without  intervening  slab  or  geometrical  optics  aporoximations.  We 
would  also  like  to  make  contact,  where  possible,  with  scattering  and 
wave  propagation  problems  in  other  areas  of  physics  in  order  to 
validate  any  proposed  solution  methods.  We  start  from  the  time 
harmonic  variation  form  of  Maxwell's  equations 

V  X  E  =  -IB  (1) 

V  X  R  =  in  =  ie(r,oj)-E  (2) 

where  e  is  the  usual  cold  plasma  dielectric  tensor.   In   this   form  of 

Maxwell's   equations  we  have   assumed   that   r  is  non-dimensional  and 

measured  in  units  of  c/w.  For  the  cold  plasma  model  e  is  hermitian  and 

with  coordinates  in  which  the  magnetic  field  is  in  the  z  direction 


-   -   (j) 


where 


v° 


-ibj    a^    0    \  (3) 
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J     Pi       c . 


and  where  to  (r)  is  the  plasma  frequency  of  the  j —  species  and  Q  (r) 
is  the  corresponding  cyclotron  frequency.  We  could  easily  include  a 
tensor  permeability  as  well  as  tensor  dielectric,  but  no  new  ideas  are 
involved  and  there  is  little  relevance  to  plasma  physics. 

Despite  the  solid  physical  basis  for  the  system  (1),(2)  it  does 
have  certain  mathematical  difficulties  that  are  most  obvious  if  we 
eliminate  B  to  obtain 

V  X  (V  X  E)  =  £(r,a))«  F,  .  (4) 

One  of  the  most  elementary  properties  of  a  set  of  differential 
equations  is  its  type:  elliptic,  hyperbolic,  parabolic,  or  mixed.  The 
type  is  usually  associated  with  solution  methods  for  the  system.  We 
may  obtain  the  type  by  looking  at  the  characteristic  polynomial  of  the 
system.  We  keep  only  the  highest  derivatives,  replace  V  by  k  and  the 
variable  E  by  6E.  We  interpret  k  as  the  normal  vector  to  some 
surface,  6 E  as  the  normal  gradient  to  the  surface  and  we  consider  the 
solvability  for  6E.^^  In  this  case  the  system  is  solvable  if  and  only 
if 
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Co  =  ^^^ 


-(k2+k2) 


KK 


k^ky 


-(kl+kh 


^x^z 


N'^z 


(5) 


S^^z 


-(k2+k2) 


is  not  zero.  Additionally,  C  is  the  symbol  of  the  operator,  and  its 
properties  describe  the  basic  nature  of  the  differential  equations,  for 
instance  the  degree  of  C  is  the  order  of  the  differential  equations. 
An  elementary  computation  shows 


C^=  0 


thus,  we  are  unable  to  Identify  the  tvpe,  or  even  the  order  of  the 
system  (1),(2).  Obviouslv,  the  system  is  highly  degenerate  and  the 
source  of  the  degeneracy  is  clear:  the  operator  Vx  applied  to  any 
gradient  vanishes.  Thus,  Vx  has  a  large  null  space,  and  this  null 
space  changes  the  character  of  the  system. 

The  introduction  of  the  usual  electromagnetic  potentials  resolves 
the  degeneracy  and  permits  direct  analysis.  We  exhibit  two  possible 
sets  of  electromagnetic  potentials  with  slightly  different  properties. 
In  general  we  may  introduce  a  vector  potential  A  and  a  scalar  potential 
(|>  so  that 


B  =  V  X  A 
E  =  -lA  -  Vij) 


(6) 
(7) 


and  (1)  is  identically  satisfied  while  (2)  becomes 


-14- 
V(V'A)  -  AA  =  e«A  -  iC'Vj    , 

the   governing  differential   equation.   As  is  well  known,  the  fields  E 
and  B  are  invariant  under  the  gauge  transformation 

A  =  A'  +  VX 
(J>  =  t  '  -  iX  . 

Thus,  we  expect  to  adjoin  one  constraint  on  A  or  ({i  in  order  to 
specify  the  potentials  uniquely. 

In   the   first   set   of  potentials,  Aj  and  >}>  ^ ,  we  select  radiation 
gauge 


iV.Aj  =  0 


(8) 


so  that  our  system  is  (8)  plus 


AA,  -  ie«V(j) ,  =  -€•  A, 


(9) 


If   we   now   repeat   the   procedure   for   the    construction    of    the 
characteristic  form  of  the  system  we  find 


k^ 
0 
0 
ik. 


0 

k2 
0 
ik. 


0 

n 


k2 


ik. 


■i(£*ls)x 


■Ke-k), 


-i(£-!s)z 


(10) 


or 
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Cj  =  -k^k*£*k  . 


(11) 


We  note  that  C^  is  not  identically  zero,  and  that  it  is  of  sixth  degree 
in  k  corresponding  to  the  reasonable  conclusion  that  (8),  (9),  or  (1), 
(2),  is  a  sixth  order  system.  We  examine  the  consequences  of  (11) 
shortly  after  we  introduce  the  second  set  of  potentials.  The  system 
(8),  (9)  is  reasonable,  but  we  see  from  (10)  that  the  system  is  not 
symmetric  or  self-adjoint. 

We  may  formulate  the  system  in  a  self-adjoint  manner  if  we   select 
a  different  gauge.   We  introduce  Aj  and  ((>  2  such  that 


iV'  (£'^'A2)  =  0  , 


(12) 


where  e"*"  Is  the  adjoint  of  e,  (e''")^^  =  e^^.   The  remaining  differential 
equation  is 


-V(V'A2)  ■*"  ^^2  "  ^1*^*2  =  "1*^2 


(13) 


For  these  potentials  the  characteristic  form  is 


k2-k2  -k^ky 


k^ky 


k2-k2 


-k  k 

K.yK.2 


i(k-e+)^   i(k-£+)y 


~  -  'x 


'^x^z 


-kyk^ 
k2-k2 

Kk-e"*"). 


-i(l*!s)x 
-i(l*b>y 

-i(£*!s)z 


(14) 


so  that 
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C2  =  -k2(lS*£'!s)(lS*£"^*lS)  •  05) 

Now,  the  form  (14)  is  clearly  self  adjoint,  and  the  characteristic 
polynomial  (15)  is  not  substantially  different  from  that  for  the  first 
set  of  potentials  (11).  We  note  that  the  self-adjoint  property  of  the 
system  does  not  depend  on  any  assumption  about  e.  The  system  (12), (13) 
is  the  set  of  Euler  equations  of  the  variational  form 

I  =  j  (dV)  {-V  X  A*«V  X  A  +  A*'£*A  +  i(j)V«(A*'e)  -  iij)  *V.  (e"*"*  A) }  .  (16) 

If  e"*"  =  e  then  the  full  variational  form  (16)  and  not  just  the 
differential  operator  is  self-adjoint.  Such  a  variational  form  should 
lead  to  good  numerical  solution  schemes  by  finite  element  or  other 
finite  difference  approximations  when  appropriate  boundary  conditions 
and  other  subsidiarv  relations  are  given.  Consistent  with  the 
variational  form  it  should  be  possible  to  specify  ^j  arid  the  tangential 
component  of  E  on  a  boundary.  Although  E  considered  as  a  solution  of 
the  basic  differential  equation  (4)  is  also  the  variational  solution  of 
(l6)  with  <\i  =<(>*=  0  when  E  is  identified  with  A,  this  latter  form  is 
degenerate  and  has  given  rise  to  serious  numerical  di  f  f  icult  ies. '^  5'  We 
expect  that  (16)  as  formulated  should  not  give  rise  to  such  problems. 

We  expect  that  a  time  harmonic  wave  propagation  problem  or  a 
scattering  problem  should  give  rise  to  an  elliptic  set  of  differential 
equations,  as  happens  in  electromagnetic  theorv  or  quantum  mechanics. 
We  shall  see  that  this  typical  situation  fails  for  the  cold  plasma 
model.   The  necessary  and  sufficient  condition  for  ellipticity  is   that 
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the   characteristic   form   (11)   or   (I5)   have  no  zeroes  for  any  real, 
non-zero,  vector  k.   The  condition  reduces  to  the  assertion  that 

C  =  k^e-k  (17) 

never  vanish  for  real,  non-zero,  k.   If  we   rotate   axes   so   that   the 
equilibrium  magnetic  field  lies  in  the  z  direction  we  find  that 


(j)     T         J  (j)     J 


or  the  necessary  and  sufficient  condition  for  ellipticity  is 

[1  -  r   (o2  /(c,2_j^2)j.[i  _  Q-^2  j;^2J  >  0  .         (18) 
(j)     J  .    1    ^ 

When  the  condition  (18)  fails  the  svstem  (8),  (9)  has  two  real 
(hyperbolic)  characteristics  and  the  system  (12), (13)  has  four  real 
characteristics.  For  frequencies  much  less  than  the  plasma  frequency 
of  the  medium,  I / i)   w^  ,  the  condition  for  ellipticity  is 


1  <  I   <.2_/(co2-n2)  .  (19) 


More  generally,  the  first  factor  in  (18)  changes  sign  when   inequality 

in   (19)   is   replaced   by  equality,  the  condition  for  what  are  usually 

called  the  lower  and  upper  hybrid   resonances.    Thus,   the   lower   and 

upper  hybrid   resonances   may   be   interpreted   as  the  occurrence  of  a 

transition  from  an  elliptic  to  a  hyperbolic  system.   Another  transition 

occurs   at   the   cyclotron  resonances  where  iji      =   H^  .   It  is  clear  that 

1 

the  cold  plasma  model  fails  near  these  values  as  o   becomes   unbounded. 
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Flnally  the  second  factor  In  (l8  )  changes  sign  at  the  plasma  resonance 
<^  =1(1)  '^t,'  The  meaning  of  the  terra  resonance  Is  quite  different 
when  we  refer  to  cyclotron  resonances  at  which  the  cold  plasma  model 
fails  or  to  hybrid  and  plasma  resonances  where  the  cold  plasma  model 
appears  to  continue  to  be  satisfactory.  In  the  latter  cases  the 
infinite  homogeneous  dispersion  relation,  slab  models,  or  geometrical 
optics  indicate  the  presence  of  "resonance"  phenomena.  It  is  our 
purpose  to  explore  further  in  subsequent  sections  the  cyclotron  and 
hybrid  resonances  and  their  significance  and  relation  to  possible 
energy  absorption  in  a  plasma. 

\  transition  from  ellipticity  to  hvperbolici ty  is  unusual  in 
mathematical  physics,  but  it  does  occur  in  fluid  dynamics  and  in  plasma 
physics,  for  instance.  The  existence  of  such  a  transition  may  give 
rise  to  singular  differential  equations  and  requires  that  we  decide  how 
to  select  solutions  from  a  potentially  singular  set.  Steady, 
compressible  fluid  flow  in  the  transonic  regime  generates  a  transition 
from  elliptic  to  hyperbolic  flow  across  the  sonic  surface.'^  The  flows 
chosen  are  smooth  across  the  sonic  surface  and  satisfy  a  causality 
condition.  For  transonic  flows  the  causality  condition  is  equivalent 
to  the  specification  of  data  on  the  real,  forward-facing 
characteristics,  no  data  on  the  backward-facing  characteristics,  and 
typical  data  in  the  elliptic  part  of  the  region.  In  the  guiding  center 
plasma  of  guiding  center  fluid  there  is  also  a  transition  from  an 
elliptic  system  to  a  hyperbolic  system.^  ■'  In  this  case,  the  transition 
is  indicative  of  instabilities  not  resolved  by  the  models.  Here,  a 
sonic  transition  indicates  the  failure  of  a  model  and  an  incomplete 
physical  representation  of  the  system. 
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Yet  another  method  of  selecting  solutions  of  wave  propagation 
problems  is  to  require  that  the  solution  be  causal  in  another  sense. 
Causality  requires  that  the  solution  [E(x,w )  ,B(x,aj )]  possess  an 
analytic  continuation  in  the  complex  (J  plane,  Im  w  <  0  and  this 
continuation  must  be  bounded  as  Im  o)  ->•  -» .  Causality  in  this  sense 
follows  from  the  observation  that  the  time  dependent  problem  with 
solution  [E(x, t ) ,B(x, t) ]  must  possess  a  well-behaved  Laplace  transform 
with  Laplace  transform  variable  p  =  iw ,  and  chat  the  time  harmonic 
dependence  solution  is  the  limit  of  the  Laplace  transform  solution  as 
Re  p  =  -Im  (jJ  tends  to  zero.  The  time  dependent  solution  Is  a 
manifestly  well-behaved  system  If  one  recalls  that  the  cold  plasma 
conductivity  comes  from  the  linearized  momentum  equation  for  each 
species 


'j%(E)  ^=  Ml5  +  -i  '^  5^"^J 


and  the  current  which  Is  to  be  inserted  In  Max'^ell's  equations 


J  =  r  e.no(r)u. 
(j) 


This  system  and  Maxwell's  equations  are  e xno ne n t i a  1 1 y  bounded, 
so  that  a  Laplace  transform  must  exist.  With  this  Interpretation 
E(x,a))  and  B(x,to)  may  be  singular  In  r  when  w  is  real.  This  approach 
has  been  applied  in  many  cases  to  plasma  wave  propagation  problems.^ 
In  particular  with  another  author  we  have  applied  It  to  lower  hybrid 
propagation  for  a  perpendicularly  stratified  medium,-  similar  to  many 
older  analogous  treatments.    We   believe   that   the   problem  of   wave 
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propagation  in  a  cold  plasma  employs  each  of  these  methods  of  selecting 
a   physically   relevant   solution  and  has   a   few  additional   unique 
questions  yet  to  be  resolved. 
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III.   WAVE  PROPAGATION  IN  AXISYMMETRIC  GEOMETRY 

In  the  previous  section  we  studied  the  cold  plasma  model  in  the 
general  case.  We  showed  that  the  wave  propagation  problem  reduced  to  a 
sixth  order  set  of  differential  equations  for  the  electromagnetic 
potentials  A  and  (fi  .  For  the  remainder  of  the  paper  we  turn  to  the 
experimentally  significant,  but  special,  case  of  wave  propagation  for  a 
plasma  equilibrium  which  is  axisymmetric.  We  show  that  with  a 
different  choice  of  potentials  the  governing  system  of  equations  is  of 
fourth  order,  and  we  show  that  these  equations  are  also  derivable  from 
a  variational  principle  different  from  (17)  when  e  is  hermitlan.  We 
examine  the  type  of  the  system  of  equations  and  compare  the  elliptic  to 
hyperbolic  transition  In  this  case  to  the  previous  case.  For 
application  In  subsequent  sections,  we  also  present  simple  explicit 
representations  of  e  and  various  pertinent  expressions  for  an 
axisymmetric  equilibrium.  We  employ  the  usual  cylindrical  coordinates 
(r,9,z),  and  orthonormal  unit  vectors  (i",9,z).  When  e  Is  referred  to 
the  basis  set  (r,e,z)  then  by  axlsymmetry  the  elements  e^^^z^^e^.^, 
etc.,  are  all  functions  of  r  and  z  alone.  We  may  then  assume  that  the 
time  varying  electromagnetic  fields  have  the  9  dependence  e^**  ,  and  we 
treat  separately  the  case  m  =  0  and  m  *  0. 

A.  The  Case  m  =  0 

If  we  take  F.q  and  Bg  as  the  two  basic  functions,  or  potentials, 
then  the  i"  and  z  components  of  (1)  and  (2)  yield 

Bj.  =  -IEq.z  (20a) 

B^  =  l(rEe),^/r  (20b) 


We 
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^r   =    ^^9.''  (21a) 

D^   =   -i(rBe),z/r    .  (21b) 

may   solve    (21)    for  E^.   and  E^    in   terras    of   Eq    and   Rq    to    find 

AE^  =   i[erz(rBe),r/r  +  t  ^^   Rq.z]    +   (erzE^^e    "  ^  vQ^  zz^^^Q  (22a) 

^^z   =  -i[trr^rnQ),^/v   +  e^r   "e.^l    +   (^zr^rt    "  ^zQ^rr^^e     .  (22b) 


where 


^     =   ^rr^^^    -    ^r.e.r     .  (23) 


rr    zz  rz    zr    • 


The   9    components    of    (1),(2)    are    as    yet    unsatisfied    and    they    implv 


"^r.z   -   ^'z.r   =   -^%  (2^a) 


^r,z   -   ^z.r   =    il^9    =    i(^9r'^T   +  ^99  1^^    +  ^gz^^^z^    •         (2^^) 


Thus,  (24)  together  with  the  definitions  (20), (22)  and  (23)  constitutes 
two  second  order  equations  for  Eg  and  Rg  .  It  is  easv  to  construct  the 
characteristic  form  C  of  the  system  and  we  find 


C  =  k2k.e'k  (25) 


where  k  =  (kj-jO.k^).  Clearly  (25)  is  a  reduction  of  the  previously 
given  forms  (11)  or  (15)  to  the  case  with  kg  =  n  corresponding  to  the 
hypothesis  of  no  9  variation  of  the  several  functions.  The  forms  that 
characterize   ellipticitv,   (18)  and  (19),  no  longer  apolv  as  v;e 
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can  no  longer  rotate  vectors  at  will  as  was  required   to  obtain   these 
forms.    Rather   than  write   out   at  this  point  the  explicit  equations 
corresponding  to  (24)  with  the  substitutions  (20),  (22)   and   (23),   we 
give  a  variational  form  whose  Euler  equations  generate  this  sytem, 

I  =  i  rdrdz  {  [V(rEe)-7(rEe)]/r2  +  [V  (  rBj)-e.V  (rBg  )  ]  /  (r^A  ) 

-  iEj[(rB9),^(erzeer  "  ^Sz^rr)/-^  +  ^  .z^^Br^zz  "  ^zr^ez^l^^ 

-  bJbq  -  EjEgfdet  (e)l/A}  (26) 

In  the  form  (26)  Vj)  =  r^ ,  j.  +  z^,^.  If  e  is  hermitian  then  (26)  is 
self-adjoint.  Such  variational  forms  frequently  form  the  basis  for 
efficient  computational  algorithms. 

B.  The  Case  m  *  0 

The   general   case   is   much  more  difficult  to  treat.   We  start  by 
considering  the  prototype  equation 

VxL  =  iM  .  (27) 

Now,  L  and  M  satisfying  (27)  exist  only  if 

V-M  =  0  .  (28) 

Hence,  we  invert  the  role  of  the  equations  (27), (28)  and  we  assume  that 
(28)  holds  and  ask  for  solutions  to  (27),  It  is  easy  to  see  that  if  M 
satisfies  (28)  then  a  particular  solution  of  (27)  is 


-24- 


L'  =  (-rM^/m.O.rM^/ra) 


so  that  the  most  general  solution  of  (27),  given  (28)  is 


L  =  (-rM^/m.O.rM^/m)  +  Vx  . 


(29) 


Hence,  the  equation  (27)  is  replaced  by  (28)  and  the  representation 
(29).  We  now  apply  this  approach  to  our  original  system  (1),(2)  and  we 
have  the  two  differential  equations 


V.  B  =  0 
V.  D  =  0 


(30) 
(31) 


and  the  representations 


E  =  (rB2/m,n,-rBp/m)  +  VX 
B  =  (-rD2/m,0,rDj./m)  +  Vu 


(32) 
(33) 


The   9    components    of    (32)    and    (33)    are 


E0    =    i  mX  /  r 
Bq    =    imu/r    , 


(34) 
(35) 


so  that  the  potentials  X  and  u  are   essentially  Eg  and   Bg .    Thus, 

whatever  value  m  has,  the  basic  potentials  are  Eg  and  Bg . 

The   r   and   z   components  of  (32)  and  (33)  are  equivalent  to  the 
system  for  D^,  D^,  Bj.  and  B^ 
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D,   -   (r/m)(e^^B^   -  e^.B^    =    ie^Q"^/'^  +^zr^.r  ^^zz^'z      ^^6) 
Bj.  +   (r/ra)D2   =  ^  'r 
Bz   -    (r/m)Dr   =  u ,^ 

We  may   easily   solve    (36)    for   D^.,    D^,    B^,    and   B^    in   terras   of  X    and  u      Co 
find 

^D^   =    i(m/r)[€^    +   (r/m)2(e^Qe^2   -  e  ^c  ^^)]\    +    [e  ^j.   -    (r/m)2Al    [\  ,^  +   (r/m)v,^] 
+  ^rzf^.z   -    (i-/n')u,rl 

TfD^    =    i(m/r)[e2e    -^   ^"^Z™)^^^  r«^  zr  "  ^zS^rr^^'^    "^  ^zr^-*^  'r  ^   (r/"')^'  .zl 
+    [e^z  -    (r/m)2A][X,^  -    (r/m)u  ,^] 

"SBj.   =   -ile^e    +   (r/m)2(ej^e2r  ~  ^zB^rr^'^    "   (r/m)£  ^^^[X  ,j.  +   (r/m)^,^] 
-   (r/m)[e^^   -    (r/m)2A]X,^   +    [1    -    (r/m)2e^Ju,^  (37) 

TTB^   =    i[e^    +   (r/ra)2(e^Qe^^  -  e^e^^)]X    +   (r/m)[e^j.  -    (r/m)-A]X,^ 
+   (r/m)   Ej-zt^'z   "    (r/m)y,rl    +    [1    -    ( r/m)  ^e  ^.^.Jp  ,^    , 


where 


T  =    1    -    (r/ra)'^(e^^4€^^)   +   (r/m)''A    .  (38) 
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To  complete  the  representation  of  the  fields  in  terms  of  the  potentials 
we  need  only  give  Dg  ,  but  Dq  =  £9^.^^  +  GggEg  +  ^62^2  so  that  with  (32), 
(34)  and  (37)  we  find 

TTDq  =  i(m/r)(eQQ  +  ( r/m)- [^qq  (e  ^.^.-fc  ^2)  -  eg^e^  -  eg^e^]  +  (r/m)^  det  e }  X 
+  [Eg^  +  (r/m)2(eQ^e^^-eQ^e^^)l  [X,^  +  (r/m)M,J 
+  [eg^  +  (r/m)2(eg^E^^-eg^e^j.)]  [X,^  -  (r/m)u,^I  .   (39) 

We  now  see  that  n  and  B  are  represented  by  the  potentials  by  (35), 
(37)  and  (39),  so  that  the  governing  differential  equations  are  (30) 
and  (31).  It  is  easy  to  show  that  the  system  consists  of  two  second 
order  partial  differential  equations  for  X  and  y  and  that  the 
characteristics  of  the  system  are  given  by  the  same  form  as  the  case 
m  =  0,  (25).  It  is  clear  that  our  equations  are  Euler  equations  of  the 
variational  form 

I  =  -  J  rdrdz(X*V.  n  +  u*V*H)  , 

or 

I  =  i  rdrdztDj.(rX*),^/r  -  DQX*(im/r)  +  ^z^*,^ 

+  B^(ru*),j./r  -  Reu*(im/r)  +  "^^v^r^]    .   (40) 

It   is   tedious   but   straightforward   to   show   that   (40),   with   the 
substitutions  (35),  (37),  and  (39)  is  self-adjoint  if  e  is  hermitian. 
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C.  Explicit  Representations  of  e  and  Associated  Functions 
Before  we  turn  to  applications  of  the  material  just  presented,   we 

give   the   explicit  structure  of  the  e  matrix  and  associated  functions. 

We  may  take  the  equilibrium  magnetic  field  to  have  the  form 


Bq  =  Bgie  cos  a  +  sin  a ( r  sin  B  +  z  cos  6)]  ,     (41) 


where  a  and  3   are   functions   of   r  and   z   alone.    \     right-handed 
orthonorraal  set  of  vectors  is 

Vj  =  r  cos  &    -   z   sin  6 

Vo  =  9  sin  a   -   cos  a ( r  sin  6  +  z  cos  &  ) 

Vt  =  0  cos  a  +  sin  a(r  sin  B  +  z  cos  6 ) 

in  terms  of  which 


e  =  I  -  (v,v,+V2V2)  1      X./(l-Y^)  +  i(vjV2-V9Vi)   T  X.Y./(1-Y^)  -  ^3^3  "[     X. 

^     ~  (1)  ^   (i)  (j) 


where 


X.  =  0)2  f^2  (42) 

J    Pj 

Yj  =  n^_/to  .  (43) 


We  have  thus  obtained  the  components  of  the  e  tensor; 
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where 


rr 


'96 


zz 


re 


a  -  b  sin  a  sin 


a  -  b  cos"  a 


7  J 

a  -  b  sin  a  cos"^ 


eg _  =  -b  sin  a  cos  a  sin 


-  1  c  sin  a  cos 


^9z  ~  ^z9 


-b  sin  a  cos  a  cos  6  -  1  c  sin  a  sin 


rz 


G^j.  =  -b  sin  a  sin  g  cos 


+  i  c  cos  a 


(44) 


a  =  1  +  )_   X./(Y^-1) 
(j) 


b  = 


r  X.[l  +  1/(Y^-1) 
(j) 

a  +  )."  X.  -  1 
(1) 


c  =  I      X^Y^/(Y3-1)  . 
(j) 


(45> 


Other  functions  of  e  which  we  need  for  the  various  representations 


are 


^  =  ^rr^zz  -  ^rz^zr  =  a(a  -  b  sin^  a)  -  c'^   cos^  a 


^1  "  ^(^ze^rz~^re^zz^  "  c(b-a)  sin  a  cos  B  -  i(c--ab)  sin  a  cos  a  sin  6 

f9  =  i{t^z^^-z^QZ^^)    =    -c(b-a)  sin  a  sin  B  -  i(c--ab)  sin  a  cos  a  cos 

h  =  ^re^Or  +  ^ze^ez  -  ^ee^^rr+^zz)  =  2a(b-a)  +  (c-ab)  sin^  a 

^4  =  ^rr  +  ^zz  =  2a  -  b  sin^  a 


(46) 


and  det  e  =  (c^-a2)(b-a) 
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The   necessary  and   sufficient  condition  for  elllpticity  is  that  k«e«k 
not  change  sign  or  that  the  form 

k.2(a  -  b  sin^  a  sin^  3)  -  2kj.k2  b  sin^  a  sin  6  cos  B  +  k^(a  -  b  sin^  a    cos^  B) 
be  definite.   Equivalently  the  system  is  elliptic  if  and  only  if 


a(a  -  b  sin^  ct)  =  a[a  cos^  a   -   (   I     X.  -  1)  sin^  a]  >  0   (47) 

(i) 


When  (47)  fails  the  system  has  two  real  (hyperbolic)  characteristics. 
The  condition  for  ellipticity  in  the  three  dimensional  case,  (18),  is 

a(a-b)  >  0  .  (48) 

If  (48)  holds  then  it  is  direct  to  show  that  (47)  holds,  or  if  the 
system  is  elliptic  in  the  three  dimensional  sense,  then  it  is  elliptic 
in  the  two-dimensional  sense.  If  (48)  fails  then  (47),  may  hold  or  it 
may  fail,  that  is,  two  dimensional  sense  ellipticitv  does  not  imply 
three  dimensional  sense  ellipticity. 

As  is  well-known,  the  quantities  a,b,  and  c  are  grossly  disparate 
in  magnitude  for  applied  frequencies  less  than  or  on  the  order  of  the 
ion  plasma  frequency.  In  the  next  two  sections  we  examine  wave 
propagation  at  frequencies  on  the  order  of  the  Ion  cyclotron  frequency 
and  on  the  order  of  the  Ion  plasma  frequency  and  lower  hybrid 
frequency.  We  look  to  see  if  the  system  is  slgnl f leant Iv  simpllfed  in 
the  frequency  range  In  question.  We  also  examine  the  implications  of 
resonances  —  or  elliptic  to  hyperbolic  transitions,  and  other 
potential  singularities  of  the  equations.  The  basic  small  parameter 
used  scale  other  quantities  is  the  mass  ratio  n. 
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n  =  m_/ra+  (49) 

and  provided  we  assume,  as  is  typical  of  many  laboratory  plasmas,  that 


^pe/f^ce  =  0(1)  (50) 


then   the   other   non-dimensional    frequency   ratios   are 

'^pe/'^pl  =    l/'^^ 

■^pe/^ci  =  0(l/n) 

"ce^^pl  =   0(1 /"^n") 

"ce/^cl  =    1/^ 

tOpi/n^i  =  o(i/A)  .  (51) 

We  have  now  laid  out  the  basic  material  we  need  for  applications. 
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IV .   WAVES  AT  ION  CYCLOTRON  FREQUENCIES 

In  this  section  we  examine  the  nature  of  the  wave  propagation 
problem  at  frequencies  on  the  order  of  the  ion  cyclotron  frequency.  We 
examine  the  resonances,  elliptic-hyperbolic  transition,  other 
singularities  of  the  system,  possible  simplifications  of  the  system, 
and  necessary  rescaling  of  the  variables.   When  we  assume  that 

n.  /o)  =  0(1)  (52) 

■^i 

it  follows  easily  from  the  definitions  and  (51)  that 

X^  =  0(l/n) 

Y^  =  0(1)  (53) 

Xg  =  0(l/n2) 

Yg  =  0(l/n) 

so  that 

a  =  0(l/n) 

b  =  0(l/n2)  (5A) 

c  =  0(l/n) 

The  problems  in  a  mirror  machine,  cos  a  =  0,  and  a  Tokamak,  cos  a  ^  1, 
are  very  different  in  structure.  Another  potential  complication  in  a 
Tokamak  concerns  the  ratio  of  the  poloidal  magnetic  field  to  the  total 
magnetic  field,  which  is  sin  a  and  is  small.  In  the  frequency  range  of 
investigation  even  sin^  a,  which  is  on  the  order  of  a  few  precent  or 
more,  is  much  larger  than  n .  Thus,  we  do  not  treat  sin  a  or  sln^  a  as 
a  small  parameter  to  be  scaled  in  n.   In   the   next   section   on   lower 
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hybrid   frequencies   we   will   have   to   scale  sin  a  with  n  in  order  to 

maintain  physical  relevance. 

We  first  examine  the  elliptic-hyperbolic  transition  and   resonance 

in   the   given   frequency  range  (52).   For  a  mirror  machine  cos^  a  =  0, 

a"d  1(^^)   Xj  >>  1,  so  that  the  system  is  elliptic,  see  (47),   for   a  <  0 

or  w^  >  Q^f  .    Thus,   above  cyclotron  frequency  the  system  is  elliptic 
i 

and  below  that  frequency  it  is  hyperbolic.  Here,  the  elliptic 
hyperbolic  transition  and  cyclotron  resonance  coincide.  Por  a  Tokamak 
the  transition  is  different.  I'Jhen  a  <  0,  (47)  clearly  holds  so  that 
the  system  is  elliptic.   When  a  >  0  the  inequality  (47)  requires 


7 

a  c 


os^  a  >  [  I     X.-Dsin-^  a  , 
(j) 


or  to  leading  order  in  n 


Yt  -  1  <  (m^/m^)  cot-  a 


The  elliptic-hyperbolic  transition  occurs  at 


J^^  /w  =  1  +  (l/2)(mg/m^)  cot-  a  + 


(55) 


a  frequency  slightly  different  from  the  cyclotron  frequency.  In  both 
cases  the  system  is  elliptic  for  high  frequencies  and  hyperbolic  at 
lower  frequencies. 

We  examined  the  situation  for  a  mirror  machine  near  resonance  In 
another  paper  (1).  In  this  case  the  resonance,  elliptic-hyperbolic 
transition,  and  breakdown  of  the  cold  plasma  model  all  coincide.  There 
we  started  with  the  Vlasov  equation  and  used  a  boundary  laver  expansion 
of  the  system  near  resonance.  We  showed  that  Che  cold  plasma  model 
could  still  be  used  provided  we   Imposed   appropriate   jump   conditions 
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across   the   resonance.    On  either  side  of  the  resonance  the  electric 
field  had  the  structure 

E  =  E^°^(y)  +  X  log  Ixl  E^^hy)   +   x  E^^Hy)   +  ...  ,    (56) 

where  x  is  a  coordinate  orthogonal  to  the  resonant  surface  and  y  is  a 
coordinate  in  the  surface.  The  electric  field  also  satisfied  the  jump 
conditions  across  the  resonance  surface 

[E]  =  0 

[(r\*V)E]  =  S«E      ,  (57) 

where  S  is  a  specified  matrix.  More  importantly,  we  could  obtain  the 
relations  (56), (57)  solely  within  the  cold  plasma  model  by  the 
imposition  of  one  form  of  causality  condition  on  the  solution.  We 
recall  that  for  causality  to  hold  in  this  case  the  solution  is  obtained 
by  the  limiting  procedure  Ira  u  +  D,  Im  cj  <  0. 

In  a  Tokaraak  there  is  an  elliptic  hyperbolic  transition  (55)  and 
the  possibility  of  a  nearby  cyclotron  resonance.  We  note  that  for  a 
very  low  beta  Tokamak  cot^  a  can  be  quite  large,  so  that  the  transition 
may  be  well  away  from  cyclotron  resonance.  Further,  if  one  uses  the 
variational  forms  for  the  wave  propagation  equations  together  with  the 
many  defintions  and  explicit  forms  (45)  and  (46)  it  is  direct  to  show 
that  the  forms  (26)  and  (40)  are  well-behaved  as  w  ->  ~^c'  Thus,  the 
only  singularity  of  the  system  is  the  elliptic-hyperbolic  transition, 
which  typically  is  very  close  to  cyclotron  resonance.  Even  though  the 
only  singularity  occurs  when  (55)  is  satisfied  we  know  that  the  cold 
plasma  model  fails  for  cj^  =  q2  ^  Except  In  a  very  low  beta  Tokamak  we 
expect  the  transition  (55)  to  occur  in  a  boundary  layer  region  in  which 
the  Vlasov  equation  must  be  employed,  exactly  as  we  did   for   a   mirror 
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machine.    We   believe   that  if  such  an  analysis  is  carried  out  answers 
very  similar  to  (57)  would  be  found.   We  do  not   present   such   results 
here. 

There  is  one  other  pair  of  possible  singularities  in  the 
differential  equations.  If  we  examine  (22)  for  the  case  m  =  0  and  (37) 
for  the  case  ra  *  0,  it  is  clear  the  the  representations  are  singular 
for  A  =  0  when  m  =  0  and  TT  =  0  for  m  *  0,  Both  singularities  occur 
when  one  solves  for  the  fields  in  terras  of  the  potentials.  We  believe 
that  the  singularities  are  only  apparent  and  do  not  show  up  in  the 
solutions  in  a  significant  manner.  In  the  full,  three  dimensional 
representation  of  the  solutions  neither  A  =  0  or T  =  0  gave  rise  to  any 
singular  behavior  of  the  solutions.  It  was  only  with  the  particular 
representations  in  terms  of  potentials  that  any  singularities  occur. 
Thus,  we  expect  that  the  solutions  for  the  potentials  will  be  such  that 
E  and  B  are  well-behaved  near  A  =  0  or  "S"  =  0.  We  expect  the  solutions 
for  Eg  and  Bg  to  be  regular  near  those  surfaces  and  such  that  E  and  B 
defined  by  (20), (22)  or  (32), (33),  and  (37)  may  be  defined  in  the  limit 
as  A  or  "S"  tends  to  zero.  Further  examination  of  the  effects  of  these 
apparent  singularities  may  be  carried  out  by  the  methods  of  the 
following  section. 

We   may   finally  consider  whether  the  scalings  (52),  (53)  and  (54) 

yield  a  significantly  simpler  set  of  equations  than  (22),  (23),  (24)  or 

(30),   (31),   (35),   (37),   (38),   (39).    We   recall  that  we  assumed  a 

distance   scale   length   c/w  =  0(c/J^   ).    For   the   waves    at    these 

i 

frequencies   the   correct   scale   length  is  c/w   ^  a  much  smaller  unit. 
Thus  i f  we  introduce 


r  =  r  f2^  /w„ 
c^   Pi 

i  =  z  f^c./^p   .  ^58) 


-35- 
then  we  expect  derivatives  with  respect  to  r  and  z  to  be  0(1).    It   is 
simpler,  and  consistent  with  (58)  to  introduce 

r'  =  /^  r 

z'  =  /n"  z  (59) 

and  again  derivatives  with  respect  to  r'  and  z'  are  0(1)  in  ri  .  It  is 
now  tedious,  but  straightforward  to  introduce  the  scalings  (52),  (53), 
(54)  and  (59)  into  the  governing  differential  equations  and  provided  we 
take 

Ee  =0(1) 

Bq  =  0(n"^/2  (60) 

when  m  =  0  and 

X    =  0(1) 

p  =  0(ri"^/2)  (f^l) 

when  m  *  0,  we  see  that  all  the  coefficients  are  essentially  of  the  same 
order   of   magnitude.   Thus,   the   scalings   leave   two   second  order 
differential  equations  with  the  same  structure.   The   coefficients   may 
simplify  somewhat,  but  no  significant  changes  occur. 

We  have  carried  out  our  study  of  ion  cyclotron  waves  without 
explicitly  imposing  the  condition  that  the  component  of  E  parallel  to 
the  equilibrium  magnetic  field  must  be  small.  It  is  Interesting  to  see 
how  this  property  of  the  solution  appears  in  our  treatment.  \  direct 
evaluation  of  E||  from  (22)  or  (32)  and  (37)  in  terms  of  the  appropriate 
potentials  shows  that  Ey  is  small,  whatever  the  potentials  are, 
provided  (60)  or  (61)  hold.   Hence,  we  need  not  consider  E|,  further. 
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The  scalings  (52),  (53),  (5A)  together  with  (59)  move  the 
singularities  A  =  0  or  T  =  0  even  further  out  of  the  range  of  practical 
interest,  so  there  Is  little  reason  to  study  them  in  detal  except  very 
near  ion  cyclotron  resonance.  Near  cyclotron  resonance  more  delicate 
analysis  based  on  the  Vlasov  equation  is  needed. 

Finally  we  observe  that  the  numerical  calculation  of  wave 
propagation  in  a  Tokamak  or  mirror  machine  remains  a  difficult  problem 
unlike  other  wave  propagation  problems.  Kven  after  the  introduction  of 
the  scalings  (60)  and  (61)  and  the  jump  conditions  (57)  appropriate  to 
a  mirror  machine  or  Tokamak,  we  are  left  with  differential  equations 
which  are  partly  elliptic  and  in  other  regimes  partly  hyperbolic. 
Although  there  is  some  experience  with  such  problems  in  transonic  flow, 
there  is  little  known  about  the  actual  problems  we  have.  Presumably 
the  variational  forms  (26)  and  (40)  would  be  the  basis  of  such 
calculations . 
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V.  WAVES  AT  LOWER  HYBRID  FREQUENCIES 

In  this   section  we   turn   to  wave   propagation   problems   at 

« 

frequencies  on  the  order  of  the  lower  hybrid  —  or  ion  plasma  — 
frequency.  We  obtain  the  scaled  equations  and  we  study  the  nature  of 
the  solutions  near  the  elliptic-hyperbolic  transition,  corresponding  to 
the  lower  hybrid  resonance.  We  concentrate  the  treatment  on  the  case 
of  ra  =  0  modes  in  a  Tokaraak,  although  the  methods  apply  quite 
generally.  The  formulas  are  much  simpler  in  this  case.  We  are 
interested  in  the  properties  of  the  solution  when  we  assume 


(^p  /o)  =  0(1)  (62) 


so  that 


Xj  =  0(1) 

Y^  =  0(/n')  (63) 


Xg  =  0(l/n) 


Yg  =  0(l//n) 


and 

a  =  0(1) 

b  =  0(l/n)  (64) 

c  =  0(l//n") 

The  condition  for  elliptlcity  of  the  system,  (47)  reduces  to  a  <  0,  or 


1  +  <^S  1^1      <   w2  /u,2  ,  (65) 

Pg   Cg      Pj^ 
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whlch  we  may  write  as 


l/o)^  >  1/u)^  +  l/|n^  n„  I  .  (66) 

Pi         ^1  '=8 


Equality  in  (65)  or  (66)  is  just  Che  condition  for  the  lower  hybrid 
resonance,  so  that  here  the  elliptic-hyperbolic  transition  is  just  the 
lower  hybrid   resonance.   For  many  equilibria  Q     f2    changes  much  less 

Cj  Cg 

than  0)^  ^  go  that  If  the  transition  occurs,  the  system  is  hyperbolic  in 
the  lower  density,  outer  region  and  elliptic  in  the  higher  density, 
inner  region.  This  arrangement  is  not  at  all  typical  for  wave 
propagation  problems  and  may  well  affect  accessibility  and  antenna 
design. 

Unlike  the  frequency  range  near  the  ion  cyclotron  resonance,  where 
one  might  argue  that  the  cold  plasma  approximation  is  inappropriate, 
near  the  lower  hybrid  resonance  there  is  no  obvious  argument  that  the 
cold  plasma  model  is  poor.  Thus,  we  believe  that  we  should  be  able  to 
study  waves  in  this  frequency  domain  within  the  cold  plasma  model 
merely  by  the  imposition  of  an  appropriate  causality  condition.  For 
simplicity  we  treat  only  the  case  m  =  0,  and  ultimately  only  for  a 
Tokamak.  There  are  several  complications  that  we  must  address,  but 
first  we  give  the  explicit  equations  derivable  from  (26) 

[(rEQ,^/r],r  +  E9,zz  +  Ee  (det  e)/A  =  [f5(rBg),^  -  f*  %  ,z] /^      (67) 


[ic^^/^)(rnQ),^/r  +   {z^^|^)%,z\,^   +  [(e^2/A)(rBe  ),r/'^ 

+  (e^^/A)Be,z],^  +  Be 
=  [-f2(rE9),r/r  +  fj  F.q.zI/A  .  (68) 
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If  we  now  introduce  a  vector  in  the  direction  of  the  poloidal  magnetic 
field 

u  =  sin  3  r  +  cos  B  z  (69) 

and  an  orthogonal  vector  in  the  poloidal  plane 

V  =  +COS  6  r  -  sin  6  z  (70) 

so      that      v,e,u   form  a    right-handed   orthogonal    set,    then  we   may   rewrite 
the   system  as 


[(rEQ.^)/r],p   +   Eq.zz   +   Eg  (det   e)/A    =   {-c(b-a)    sin  a    [u'V(rRQ)]/i 
+   i(c2-ab)    sin  a    cos  a    [v«V  (rBg  )  ] /r} /A  (71) 


rV.{  [a/(r2A)]V(rBe)}    -   rV •  {  [(b   sin^  a)/(r2  A)]    u(u-V)(rBe)} 
-   i6«V(rBe)  X   V(  (c   cos  a)/(rA)}    +   Bq 

=   {r(b-a)    sin  a  [u«V  (rEg  )  ] /r  +    l(c^-ab)    sin  a    cos  ot  [v.  V  (tEq  )  ] /r} /A 

(72) 

We  see  easily  that  the  second  order  differential  operator  in  (71)  is 
essentially  the  Laplaclan,  while  in  (72)  the  First  of  the  two  second 
order  operators  is  also  the  Laplaclan  while  the  second  of  the  two  is 
essentially  (u»7)^.  The  first  equation  is  always  elliptic  while  the 
second  equation,  since  b  >>  |a|  >  0  is  elliptic  for  a  <  0  and 
hyperbolic  otherwise. 

With  the  equations  displayed  explicitly,  we  wish  to  examine  two 
further  questions.  We  wish  to  determine  the  scaling  of  the  solutions 
in  n   that   generates  reasonable  electromagnetic  fields  and  we  wish  to 
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study  the  solutions  near  the  elliptic-hyperbolic  transition  —  or  lower 
hybrid  resonance.  We  must  distinguish  carefully  the  case  of  a  Tokaraak, 
cos  a  ^  1,  and  the  case  of  a  mirror,  sin  a  =  1.  It  is  clear  that 
different  terras  of  different  orders  of  magnitude  in  n  appear  in  the  two 
cases.  We  choose  to  examine  the  case  of  a  Tokamak,  the  mirror  could  be 
done  similarly.  Tn  the  previous  section,  when  we  treated  the  Tokamak 
we  took  sin^  a  =  B^/B  to  be  of  order  one,  but  here  the  relative 
magnitudes  of  the  terms  in  (71)  and  (72)  are  not  so  different  that  we 
can  afford  to  ignore  the  smallness  of  sin^  a.  Hence,  consistent  with 
situations  of  physical  interest  we  take 

sin  a  =  0(ri^'''^)  .  (73) 

Ue      now   turn   to  the  scaling  of  the  solutions  in  n .   It  is  easy  to  see 
that  if  we  select  a  scaling  of  the  electromagnetic  wave  with 

Ee  =0(1) 

Bq  =  0(n"^/2)  (74) 

V,  =  0(11"^/^) 


Vj^  =  0(n"l''2) 


then  we  have  a  self-consistent  system  with  essentially  all  the  terms  of 

(71),   (72)   Included  and   of   the  same  order  of  magnitude.   We  do  not 

modify  the  characteristic  equilibrium  scaling  length  from  c/aj„  . 

Pi 

It  is  possible  to  see  in  (71),  (72)  with  the  scalings  (73),  (74)  a 
significant  failure  of  conventional  geometrical  optics  previously 
Identified  in  slab  geometry . '"^^  There  is  a  term  in  (72)  proportional  to 
Y,[(c  cos  6)/(rA)]  even  though  we  have  assumed  that  the  equilibrium 
varies  on  the  distance  scale  c/u   while  the  waves  vary  on   the   faster 
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distance   scale  n~  '  c/aj„  or  n~   c/(iJ„  .   Thus,  the  condition  that  the 

Pi  Pi 

waves  vary  more  rapidly  than  the  equilibrium  is  not  adequate  to  justify 
the  use  of  geometrical  optics  on  the  original  system  (1),  (2).  The 
result  should  not  be  surprising  since  the  elements  of  a  in  the  equation 
depend  non-trivially  on  the  geometrical  optics  parameter.  Thus, 
cut-off  and  propagation  properties  of  lower  hybrid  waves  derivable  from 
the  infinite  homogeneous  medium  are  incorrect  as  they  miss  such  terms. 
We  have  already  shown  In  slab  geometry  a  major  modification  in 
propagation  characteristics  that  originates  with  such  a  term.  We 
could,  however,  use  (71),  (72)  with  the  scallngs  (73),  (74)  as  the 
basis  for  a  geometrical  optics  treatment  of  lower  hybrid  waves,  and 
this  analysis  would  contain  the  effects  of  spatial  gradients.  That  is, 
we  could  take  the  system  (71),  (72)  and  assume  that  (rEg)  and  (rBg) 
vary  rapidly  in  r  and  z  while  all  equilibrium  quantities  and  their 
derivatives  vary  of  order  one  in  r  and  z.  This  system  would  then  be  the 
basis  for  the  usual  geometrical  optics  procedure.  The  introduction  of 
the  scaling  (74)  and  coordinates  parallel  and  perpendicular  to  B  before 
the  geometrical  optics  approximation  would  simplify  the  process 
considerably. 

The  properties  we  have  just  given  applv  equally  well  for  the  modes 
with  m  *  0.  One  can  show  that  the  system  of  equations  (30),  (31)  with 
the  substitutions  (35),  (37),  (39)  also  distinguishes  parallel  and 
perpendicular  derivatives  of  the  waves  and  leads  to  scallngs  of  the 
form  (73),  (74)  and  to  the  non-trivial  appearance  of  spatial  gradients 
of  equilibrium  quantities.  We  do  not  pursue  this  question  or  wave 
propagation  in  a  mirror  further  here. 

The  last  question  we  examine  is  the  nature  of  the  solution  near 
the  elliptic-hyperbolic  transition.  We  start  from  the  system  (71), 
(72).  Since  very  little  is  known  in  general  about  systems  of  equations 
which  change  type,  we  recognize  that  many  of  our  assertions  are  in  fact 
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conjectures.  We  base  our  comments  on  the  Information  available  to  us. 
We  first  transform  the  equations  to  normal  form  for  a  second  order 
system.  We  find  that  there  are  two  fundamentally  different  situations. 
When  the  resonance  surfaces,  a  =  0,  coincide  with  flux  surfaces  the 
plasma  behaves  like  a  perpendicular  stratified  medium  and  there  is 
energy  absorption  along  the  entire  resonance  —  or  flux  —  surface. 
The  results  in  this  case  could  be  obtained  by  rather  standard 
inner-outer  boundary  layer  analysis  and  are  rather  compelling.  The 
situation  when  the  resonance  surfaces  do  not  coincide  with  flux 
surfaces  Is  quite  different.  In  the  neighborhood  of  these  points  where 
the  two  surfaces  are  not  tangent  the  normal  form  of  the  equation  (72) 
Is  essentially  the  Tricomi  equation  of  transonic  flow.^^^  In  the 
neighborhood  of  such  points  the  likely  Interpretation  of  the  condition 
of  causality  Is  that  data  Is  given  on  some  characteristics  and  not  on 
others,   but   no  singularity  of  the  solution  —  or  energy  absorption  — 

is  expected.   Recall  that  the  local  Ohmic   energy   production   rate   is 

* 

Re  (E  "J),  which  vanishes  identically  for  the  cold  plasma  conductivity. 

Only  if  E  is  singular  is  heating  possible.  Near  the  points  of  tangency 
of  the  resonance  surfaces  and  flux  surfaces  it  is  not  possible  to 
transform  the  solution  smoothly  to  the  Tricomi  form.  We  believe  that 
energy  absorption  is  possible  near  these  points,  but  we  are  unable  to 
provide  a  definite  prescription  for  its  computation  beyond  the 
conventional  procedure  of  gi^/ip.g  o)  a  small  Imaginary  part  and  computing 
the  solution.  If  we  examine  a  simple  model  problem  that  corresponds  to 
the  case  of  tangency  of  the  two  surfaces  we  can  show  that  the  family  of 
special  solutions  corresponding  to  similarity  solutions  has  no 
singularities   on   the   resonance   surface,  but  may  be  singular  at  the 

point  of  tangency  of  the  two  surfaces.   Flnallv,  an  examination  of   the 

(9) 
results  of  a  geometrical  optics  solution    of  the  system  also  indicates 

that  at  points  at  which  the  resonance  and  flux  surfaces  are  not  tangent 
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no  peculiarity  of  the  solutions  appears  while  at  points  of  tangency  a 
resonance  behavior  does  occur.  Thus,  we  believe  that  energy  absorption 
occurs  only  near  the  points  of  tangency  of  resonance  and  flux  surfaces. 
As  we  noted  previously,  of  the  two  equations  (71),  (72),  the  first 
is  always  elliptic  while  the  second  exhibits  the  change  of  type.  Thus, 
it  is  plausible,  and  greatly  simpler,  to  consider  that  (71)  determines 
rEg  when  rBg  is  known,  while  (72)  gives  rBg  with  rEg  specified.  We 
consider  rEg  known  and  concentrate  on  (72)  for  rBg  alone.  Further, 
since  we  expect  the  highest  order  derivatives  to  be  dominant,  it  should 
suffice  to  consider  the  form 

F  =  a(r,z)(u,^j.-Ki,22)  +  b  sin^  a    (sin^  a    u,^^ 

+  2  sin  a  cos  a  u,^^  +  cos"  a   u,^^)  ,  (75) 

which   gives  the  leading  order  terms  in  (72).   We  are  interested  in  the 
structure  of  (75)  near  a  =  0.   If  we  introduce  new  coordinates 

C  =  a(r,z) 

n  =  n(r,z)  ,  (76) 

where  n(r,z)    is   unspecified   then 

F   =    laS^,^.  +  b   sin^  a(sin  a   C.j.  +  cos  a   C,z)^lu  ^^ 

+  2[a(C  ,j.n  ,^-K  ,2n  ,2)   +  b   sin^  a    (sin  ct   ?,^+cos  a   C,2)(sin  a   n,r+cos  a   n,^))u,^j^ 
+   (an^,j.  +  b   sin     a    (sin  a  n,j.+cos  a   ^,z)    ^u,^^    +  lower   order   derivatives    . 

(77) 

We   may   obtain   a   normal   form   for   F   if   we   eliminate   Che   cross 
derivatives,  or  if  we  select  n  so  that 
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aVC'Vn  +  b  sin^  a (u-V?  )(u«7n  )  =  0  (78) 

We  must  distinguish  the  cases  u'V^  =  0  and  ti'VC  *  0  .  If  u»7C  =  0  then 
5  =  const,  on  a  flux  surface  and  resonance  and  flux  surfaces  coincide. 
In  this  case  we  take  n  such  that  ^^'Vn  =  0  or  the  curves  n  =  const. 
are  the  orthogonal  trajectories  to  the  curves  ?  =  const,  and  the  flux 
surfaces.   Now, 

F  =  L(C,n)(5u  ^^  +  M(?,n)u,^^  +  lower  order  derivatives     (79) 

and  L(C,n)  and  M(C ,n )  are  regular  and  well-behaved  near  ?  =  0,  the 
resonance  surface  and  M(5 ,n )  >  0.  When  L(5 ,n )  and  M(C  ,n )  vary  slowly 
with  respect  to  their  arguments  as  they  do  here  it  is  plausible  that 
the  system  could  be  approximated  by  an  ordinary  differential  equation, 
and  then  the  system  looks  like  a  slab  model  with  a  singularity  at 
C  =  0.  Slab  analysis  indicates  absorption  at  C  =  0.^^^  We  do  not 
pursue  this  case  as  it  is  very  untypical  of  most  plasmas. 

When   (u.V)5  t    0,  we  may  use  (78)  to  solve  for  n,  and  in  this  case 
the  normal  form  of  F  is 

F  =  L(?  ,ri)(M(C  ,Tl  )u,^^  +  Cu.^f^j  +  lower  order  terms  .       (80) 

We  note  that  (78)  and  (79)  are  drastically  different  equations  as  the 
potentially  vanishing  factor  of  5  multiplies  different  derivatives  in 
each  case.  Now  (80)  is  essentially  the  differential  equation  of 
transonic  flow,  or  the  Tricomi  equation.^ ^^  'extensive  analysis  has 
indicated  that  if  one  gives  data  on  characteristics  going  Into  (or  out 
of)  the  singular  surface  ^  =  0  and  data  on  the  boundary  on  Che  elliptic 
side,  then  there  is  a  unique  smooth  solution  of  the  system.   Such   data 
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corresponds   to  the  data  we  would  give  for  one  of  the  two  equations 
(71), (72)   and   it   seems   quite   likely   the   solution   possesses   no 
singularity  on  the  resonance  surface.   Causality  here  reduces  to  the 
selection  of  the  proper  sets  of  characteristics  for  data. 

We  may  consider  a  simple  model   equation   that   has  many   of   the 
properties  that  interest  us.   We  seek  solutions  of  the  equation 

(x-y2)u,^j^  +  u.yy  =  0  .  (81) 

Clearly  the  equation  Is  elliptic  for  x  >  y  and  hyperbolic  for  x  <  y^. 
the  resonance  surface  Is  x  =  y^  and  the  analog  of  the  flux  surface  is 
the  curve  x  =  const.,  on  which  y  varies.  Unfortunately  we  can  offer 
very  little  information  on  this  equation,  but  we  can  study  one  simple 
class  of  solutions,  the  similarity  solutions  for  which 

u  =  x°'f(y2/x)  (82) 


(1-W)[a(a-I)f(y)  -  2(a-l)iif,^  +  l^  f'ppl  +  ( 2f  .p+^P  f  .^y  )  ,   (83) 


where 


P  =  y2/x  .  (84) 

We  note  first  that  p  =  1,  the  resonance  curve  Is  not  a  singular  point 
of  the  equation  (83).  Thus,  we  expect  no  singularities  In  the  family 
of  solutions  (82).  If,  for  Instance,  we  sought  similarity  solutions  of 
the  wave  equation  u,^^^  =  u,^„  of  the  form  u  =  g(x/y)  then  we  would  find 
the  light  cone  x^  =  y  as  a  singularity.  Returning  to  (81)-(84)  we 
note  that  for  |p  I  >>  1  (or  |x|  small)  f  ~  u"  or  f  ~  p"^"^  so  that  u  Is 
well-behaved  for  x  small.   For  our  family  of  similarity  solutions  there 
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are  no  singularities  on  the  resonance  surface  u  =  1  or  flux  surface 
X  =  0.  While  this  simple  analysis  is  not  even  conclusive  for  the 
simple  model  equation,  it  does  support  our  view  that  the  only  remaining 
singularity  possible  is  at  the  point  where  the  resonance  surface  and 
flux  surface  are  tangent  x  =  y  =  0.  At  this  point  singular  solutions 
allow  for  the  possibility  of  an  energy  sink,  for  the  wave,  or  a  plasma 
heating  zone.  Presumably  such  heating  could  be  calculated  by  taking 
Im  to  +  O"*". 

Our  last  approach  to  this  question  is  based  on  a  geometrical 
optics  approximation  solution  to  the  system. ^"^  Numerical  calculation 
of  rays  in  geometrical  optics  in  a  non-perpendicular  stratified  medium 
shows  that  the  rays  pass  through  the  lower  hybrid  resonance  with  no 
significant  modification  and  that  the  the  ray  wave  number  k  certainly 
does  not  become  unbounded.  This  result,  found  numerically,  can  be 
demonstrated  analytically  as  well,  if  one  solves  the  ray  ordinary 
differential  equations  near  a  lower  hybrid  resonance  surface.  We  do 
not  carry  out  this  elementary  exercise  in  analysis  as  it  is,  in  fact, 
essentially  the  same  argument  as  occurs  in  going  from  (77)  to  (80). 
The  numerical  work  and  elementary  analysis  both  support  our  contention 
that  except  near  the  points  where  resonant  surfaces  and  flux  surfaces 
are  tangent  no  singularities  or  wave  heating  occur. 

The  formulation  of  a  wave  propagation  problem  in  this  frequency 
range  is  rather  complex.  If  the  wave  frequency  is  everywhere  below 
lower  hybrid  resonance  then  the  system  of  equations  is  elliptic  and  the 
problem  becomes  a  standard  scattering  problem.  If  the  wave  frequency 
is  everywhere  above  lower  hybrid  resonance  then  the  system  of  equations 
is  part  elliptic  and  part  hyperbolic.  Incoming  and  outgoing  wave 
solutions  are  defined  by  the  forward  and  backward  facing 
characteristics  so  that  it  should  be  possible  to  prescribe  incoming 
waves   at   one   part   of  the  boundary  and  allow  outgoing  waves  at  other 
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parts  of  the  boundary.  When  lower  hybrid  resonance  occurs   in  the 

plasma,   we   still  have  incoming  and  possibly  outgoing  waves  at  part  of 

the  boundary,  but  we  also  have  energy  sinks   at   those   few  points   at 

which  the  flux  surfaces  and  resonance  surfaces  are  tangent.   A  more 

thorough  analysis  of  the  solutions  near  these  points,  akin  to  the  study 

of   similarity   solutions   (81)-(84),   is   needed   to   determine   the 

appropriate  behavior  of  the  solutions  there. 

Finally  we  note  that  in  Tokamaks  with  u^   <<  |n  fi   | ,  that  is  for 

Pi       1   e 

very  high  field  or  very  low  density  systems,  the  resonance  surfaces  are 
approximately  constant  density  surfaces  and  thus  likely  to  be 
approximately  flux  surfaces.  In  this  approximately  perpendicularly 
stratified  medium  we  expect  strong  heating  all  along  the  resonance 
surface  exactly  as  we  found  in  special  cases. 
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